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ABSTRACT 

We report the first results of a program to study the internal kinematics of globular clusters in the 
outer halo of the Milky Way. Using the Keck telescope and High Resolution Echelle Spectrometer, we 
have measured precise radial velocities for 30 candidate red giants in the direction of Palomar 13, an 
object traditionally cataloged as a compact, low-luminosity globular cluster. We have combined these 
radial velocities with published proper motion membership probabilities and new CCD photometry from 
the Keck and Canada-France-Hawaii telescopes, to isolate a sample of 21 probable members. We find a 
systemic velocity of {vr)s — 24.1±0.5 km s~^ and a projected, intrinsic velocity dispersion of CTp — 2.2±0.4 
km s^^. Although small, this dispersion is nevertheless several times larger than that expected for a 
globular cluster of this luminosity and central concentration. Taken at face value, this dispersion implies 
a mass-to-light ratio of Ty — 40]ti7 based on the best-fit King-Michie model. The surface density profile 
of Palomar 13 also appears to be anomalous among Calactic globular clusters — depending upon the 
details of background subtraction and model-fitting, Palomar 13 either contains a substantial population 
of "extra-tidal" stars, or it is far more spatially extended than previously suspected. The full surface 
density profile is equally well-fit by a King-Michie model having a high concentration and large tidal 
radius, or by a NFW model. We examine — and tentatively reject — a number of possible explanations 
for the observed characteristics of Palomar 13 (e.g., velocity "jitter" among the red giants, spectroscopic 
binary stars, non-standard mass functions, modified Newtonian dynamics), and conclude that the two 
most plausible scenarios are either catastrophic heating during a recent perigalacticon passage, or the 
presence of a massive dark halo. Thus, the available evidence suggests that Palomar 13 is either a globular 
cluster which is now in the process of dissolving into the Galactic halo, or a faint, dark-matter-dominated 
stellar system. 

Subject headings: Galaxy: halo — Galaxy: kinematics and dynamics — globular clusters: individual 
(Palomar 13) — galaxies: dwarf 



1. INTRODUCTION 

Much of our knowledge of the Galactic halo is based di- 
rectly upon observations of its globular clusters, the most 
readily identifiable halo objects. Those clusters which are 
located in the outer halo of the Milky Way are especially 
important probes of the formation and evolution of the 
Galaxy, as their ages and metallicitics provide direct con- 
straints on the duration of halo formation process and on 
the time-scale for Galactic chemical enrichment, while the 

^ Based on data obtained at the W.M. Keck Observatory, which is operated as a scientific partnership among the California Institute of 
Technology, the University of California, and NASA, and was made possible by the generous financial support of the W. M. Keck Foundation. 

^Visiting Astronomer, Canada-France-Hawaii Telescope, operated by the National Research Council of Canada, the Centre National de la 
Recherche Scientifique of France, and the University of Hawaii. 
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shape and extent of Galaxy's dark halo is constrained by 
their orbital properties. 

These distant clusters share their location in the outer 
halo with a system of nine dwarf spheroidal (dSph) galax- 
ies. These dwarfs have Galactocentric distances in the 
range 24 < Rq < 250 kpc (Mateo 1998) — an interval 
which includes 16 globular clusters (Harris 1996). To some 
extent, the dSph galaxies resemble the globular clusters in 
that they always contain a population of old metal-poor 
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stars; there arc, however, a number of important distinc- 
tions between the two classes of objects. The central stellar 
densities of the dSph galaxies are much lower than those of 
typical globular clusters, and many of the dwarfs show evi- 
dence for intermediate-age populations (e.g., Mateo 1998). 
More importantly, measurements of the central velocity 
dispersions of the dSph galaxies suggest that they are em- 
bedded in extended dark halos (Aaronson 1983) which are 
presumably non-baryonic in nature (Feltzing, Gilmore & 
Wyse 1999), although alternative explanations have been 
proposed (e.g., Milgrom 1995; Klessen & Kroupa 1998). 

By contrast, color magnitude diagrams (CMDs) of 
distant halo globular clusters show no sign of distinct 
intermediate-age populations, although there is evidence 
that some of the chisters may be younger than their coun- 
terparts in the inner Galaxy (e.g., Stetson et al. 1999). 
And while there is a preponderance of faint, low-density 
globular clusters in the outer halo (see, e.g.. Figure 1 of 
van den Bergh 1994), their central stellar densities remain 
comfortably above those of even the most centrally con- 
centrated dSph galaxies. 

Furthermore, dynamical studies of nearby globular clus- 
ters have established that, unlike the dSph galaxies, they 
appear to contain little or no dark matter apart from nor- 
mal stellar remnants such as white dwarfs and neutron 
stars (e.g., Pryor & Meylan 1993). However, a direct de- 
termination of the internal mass distribution in distant 
halo clusters has proven impossible since the faintness of 
even their most evolved red giants — coupled with their 
sparsely populated red giant branches (RGBs) and low 
surface brightnesses — make the requisite spectroscopic 
observations extremely challenging. For example, 12 of 
the 16 globular clusters with Rq ^ 24 kpc have central 
surface brightnesses below /Uy ~ 19.5 mag arcsec"^; for 
an assumed mass-to-light ratio^ of Ty = 2, these clus- 
ters are expected to have central velocity dispersions of 
apfi ^ 2 km s~^. Thus, to study their internal dynamics, 
a velocity precision of e{vr) ^ 1 km s~^ for metal-poor 
stars in the range 17 ^ V ^ 20 is required. 

It is hardly surprising that these objects have been ne- 
glected in previous radial velocity surveys of Galactic glob- 
ular clusters. The lone exception is NGC 2419, for which 
Olzewski, Pryor & Schommer (1993) measured velocities 
for 12 red giants. From this small sample, they found 
Ty ~ 0.7, among the lowest measured for any globular 
cluster and contrary to the notion that this cluster is em- 
bedded in a massive dark halo. Nevertheless, this object 
is — with an absolute magnitude of My = —9.6 (Harris 
1996) among the most luminous Galactic globular clus- 
ters and atypical of the clusters which populate the outer 
halo. Clearly, a dynamical study of a more representative 
sample of halo clusters is in order. 

There are several reasons to suspect that a search for 
dark halos in outer halo clusters might prove profitable. 
First, the characteristic mass of globular clusters is sim- 
ilar to the Jeans mass at recombination, and thus has a 
special cosmological significance; indeed, it has often been 
proposed that globular clusters may contain dark halos 
{e.g., Peebles 1984; West 1993). While there is no evi- 
dence to support this conjecture, it has in practice proven 
surprisingly difficult to rule out (Heggie & Hut 1996). For 



example, as clusters evolve in the Galactic tidal field, they 
are expected to lose mass throiigh tidal stripping a 
process which confounds efforts to reconstruct the initial 
mass distribution in these clusters using present-day ob- 
servations. Second, there are fairly compelling reasons to 
believe that some Galactic dSph galaxies are dark-matter- 
dominated, and since globular clusters are simply the next 
step down in luminosity, it is therefore reasonable to ex- 
pect that the same scenario which is usually invoked to 
explain the internal dynamics of the dSph galaxies (i.e., 
gas ejection from shallow potential wells by supernova- 
driven winds, followed by adiabatic expansion; Dekel & 
Silk 1986) might produce end-products that resemble glob- 
ular clusters if the gas ejection mechanism was particu- 
larly efficient. This possibility is especially topical in light 
of the emerging evidence from numerical simulations that 
the halos of large galaxies have been assembled hierarchi- 
cally (e.g., Klypin et al. 1999; Moore et al. 1999 and 2001; 
Cote et al. 2000), although the predicted number of low- 
mass dark halos greatly exceeds the observed number of 
dwarf galaxies in such simulations. Various means of re- 
solving this important contradiction have been proposed: 
i.e., by supressing the small-scale power in the family of 
cold dark matter models (see, e.g., Dave et al. 2001), 
or by postulating the existence of a large population of 
"dark galaxies" (see, e.g., Hirashita, Takeuchi & Tamura 
1999; Bullock, Kravtsov & Weinberg 2000). There are 
also some reasons to believe that isolated, low-luminosity 
clusters may contain significant amounts of haryonic dark 
matter. N-body simulations of evolving globular clusters 
suggest that as clusters lose mass through dynamical evo- 
lution, the fraction of their mass contained in white dwarfs 
should increase steadily (Vesperini & Heggie 1997). The 
outer halo contains a disproportionately high fraction of 
low-luminosity clusters, so a search for dynamical evidence 
of extreme white dwarf populations in these clusters might 
prove worthwhile. 

In 1998, we began a program to study the inter- 
nal dynamics of distant halo globular clusters using the 
High Resolution Echelle Spectrometer (HIRES) at the 
W.M. Keck Observatory. This program, which was de- 
signed to yield the first direct measurements of the velocity 
dispersions and mass-to- light ratios for these clusters, was 
also motivated by the need for improved measurements of 
the systemic velocities of several of these clusters (i.e., the 
radial velocities of four of the 16 clusters having Galacto- 
centric distances greater than 24 kpc are either unknown, 
or have uncertainties of more than 20 km s~^). It has been 
noted on numerous occasions that a number of the Galac- 
tic dwarf galaxies and globular clusters fall along great 
circles (e.g., Lynden-Bell 1976ab; Kimkel 1979; Majewski 
1994). Such structures would have obvious implications 
for the formation of the Galactic halo, but unambiguous 
evidence of their existence requires orbital parameters for 
the putative members. Programs are now underway to 
measure the proper motions of the distant halo globular 
clusters and dwarf satellites (e.g., Dinescu et al. 2001), 
and accurate radial velocities are a prerequisite for the 
measurement of their orbits. 

In this paper, we present the first results from this pro- 
gram. Our sample consists of seven, predominantly low- 
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concentration clusters having Galactocentric distances in 
the range 25 <^ Rq ^112 kpc. An analysis of the full sam- 
ple will be presented in a separate paper; here, we present 
our findings for a single sparse object, Palomar 13. With 
an absolute magnitude of My — —3.8 mag, Palomar 13 is 
exceeded in luminosity by all but a handful of the ~ 150 
globular clusters belonging to the Milky Way, and is one 
of the objects which has been identified on previous occa- 
sions as a possible member of a great stream (e.g., Lynden- 
Bell 1976a; Lynden-Bell & Lynden-Bell 1995). Based on 
our new radial velocities and photometry, we find both the 
structural properties and the mass-to-light ratio of this ob- 
ject to be anomalous among the Galactic globular cluster 
system. 

2. OBSERVATIONS AND REDUCTIONS 

2.1. Keck and CFHT Photometry 

On the night of 10 September, 1999, we used the Low 
Resolution Imaging Spectrometer (LRIS; Oke et al. 1995) 
on the Keck II telescope to obtain a series of BVI im- 
ages of Palomar 13. An observing log for these, and other, 
observations of Palomar 13 may be found in Table |^. In 
imaging mode, LRIS has a spatial scale of 0'.'215 pixel^^ 
and a field of view of 5'.8x7'.3. For each filter, we ob- 
tained a pair of images centered on the cluster. Exposure 
times were 10 and 300s in V, 20 and 480s in B, and 6 and 
180s in /. Conditions during the night were photometric, 
and the FWHM of isolated stars within the BVI frames 
were measured to be in the range 0'.'7-0'.'8. Images were 
bias-subtracted and then flat-fielded in the IRAF^ envi- 
ronment using sky flats obtained during twilight. Instru- 
mental magnitudes for unresolved objects in the field were 
derived using the DAOPHOT II software package (Stetson 
1993), and calibrated with observations of several Landolt 
(1992) standard fields taken throughout the night. The 
calibration equations took the form 



V 

B-V 
V -I 



V + ai — biXv + ci(6 — v) 

02 - B + C2{b - v) (1) 



03 - bsXi + C3{v - i). 



Aperture corrections were determined separately for each 
filter, and a master object list was created using the pho- 
tometry files for the long and short exposures. The final 
photometric database contains 840 objects detected with 
a minimum point-source signal-to-noise ratio of S/N — 5 
in all three filters. A finding chart for Palomar 13 con- 
structed from our V images is shown in Figure |l|. All stars 
having measured radial velocities (see §2.2) are identified 
on this image. 

Images of Palomar 13 were also obtained on the night 
of 14 July, 1999 using the 3.6m Canada-France-Hawaii 
telescope (CFHT) and the CFH12K mosaic camera (Cuil- 
landre et al. 2000). A series of three 600s exposures in 
each of the V and R filters were taken with the clus- 
ter centered on chip #8 of the CCD mosaic. Each im- 
age was bias-subtracted and flat-fielded, and then shifted 
and stacked to produce a final image for each of the two 
filters. Isolated stars in the stacked images were found 
to have FWHMs of 0'.'75 in each filter. Photometry was 
then performed using DAOPHOT II, and the photome- 

*IRAF is distributed by the National Optical Astronomy Observatories, 
in Astronomy, Inc., under contract to the National Science Foundation. 



try lists were matched. Unfortunately, conditions were 
non-photometric and it was not possible to calibrate the 
photometry directly. However, instrumental (w) magni- 
tudes for stars on chip #8 were calibrated via secondary 
standards in this field selected from our Keck/LRIS pho- 
tometry. Thus, the photometric database from our CFHT 
data consists of calibrated V magnitudes and instrumen- 
tal r magnitudes for 855 objects in this CCD field, which 
measures 7.'0xl4.'0. This single field is 2.3 times larger 
than that available from LRIS, although the photometric 
catalog is limited to the range 18.5 <^V <^ 23.5 mag (as 
opposed to our LRIS point-source catalog, which spans the 
range 15.5 <:,V <^ 24.5 mag). 

By virtue of its remote location in the outer halo, Palo- 
mar 13 should, in principal, play an important role in es- 
tablishing the chronology of halo formation. Somewhat 
surprisingly, it has been the focus of only little attention 
in this regard. CMD studies are limited to photographic 
study of Ortolani, Rosino & Sandage (1985), the CCD 
study of Borissova, Markov & Spassova (1997) and the re- 
cent photographic/CCD study of Siegel et al. (2001). In 
Figure |[ we compare our LRIS magnitudes and colors for 
stars in common with each of these three studies. We find 
good agreement with the BV photometry of Siegel et al. 
(2001), but the comparison reveals some interesting dis- 
crepancies with the earlier studies. There is evidence for 
a non-linearity in the V magnitudes of stars fainter than 
V ^ 21 in the Ortolani et al. (1985) study, while the 
Borissova et al. (1997) V magnitudes show a rather large 
scatter with a mean offset of ISV = 0.22 ± 0.06 mag. Due 
to the small number of red giants in this cluster — and 
the limited depth of previous photometric studies — the 
stars in common between these studies and ours span a 
modest range in color. For both the Ortolani et al. (1985) 
and Borissova et al. (1997) datasets, we find mean color 
offsets of A(B - y) - 0.1 mag. 

2.2. HIRES Spectroscopy 

In a series of observing runs during 1998 and 1999, we 
used HIRES (Vogt et al. 1994) to acquire spectra for can- 
didate red giants in the direction of seven distant halo 
globular clusters: Palomar 3, Palomar 4, Palomar 5, Palo- 
mar 13, Palomar 14, NGC 7492 and NGC 2419. Table | 
presents an observing log for our HIRES observations of 
Palomar 13, which were obtained during the course of four 
dark runs on the Keck I telescope. 

Candidate RGB stars in Palomar 13 were selected from 
previously published CMDs and finding charts (e.g., Or- 
tolani et al. 1985 and Borissova et al. 1997). A few addi- 
tional targets were identified on the basis of instrumental 
CMDs constructed from images obtained with COSMIC 
(Kells et al. 1998) on the Palomar 5m telescope. Dur- 
ing three of the four HIRES runs, we limited the entrance 
aperture to 0'.'86x7'.'0 with the CI decker and binned the 
2048x2048 detector 1x2 (i.e., in the spatial direction) to 
reduce the read noise. For the fourth run, we used the C5 
decker (i.e., an entrance aperture of dimension l'.'15x7'.'0) 
and binned the detector 2x2. The corresponding spec- 
tral resolutions for these two instrumental configurations 
are A/AA — 45,000 and 34,000 respectively. In all cases, 

which are operated by the Association of Universities for Research 
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we used a single readout amplifier, a gain setting of 2.4 
ADU^^, the red collimator and a cross-disperser in 
first order. The angles of the grating and cross disperser 
were adjusted to give complete spectral coverage over the 
range 5055 ^ A ^ 5355 A. Thorium- Argon comparison 
lamp spectra were acquired frequently during each night 
— usually before and after each program star observation. 
In a few cases, lamp spectra were separated by 30-60 min- 
utes for short exposures of the brightest program stars. 
High-S/N spectra for 8-12 lAU radial velocity standard 
stars were obtained during each observing run. 

All spectra (i.e., program objects and standard stars) 
were reduced in an identical manner following the general 
procedures described in earlier papers (e.g.. Cote 1999; 
Cote et al. 1999). The radial velocity of each program 
object was measured by cross-correlating its spectrum 
against that of a master template created during each run 
from the observations of lAU standard stars. In order to 
minimize possible systematic effects, a master template 
for each observing run was derived from a similar, and in 
some cases identical, sample of lAU standard stars. From 
each cross-correlation function, we measured both Vr, the 
heliocentric radial velocity, and Rtd, the Tonry & Davis 
(1979) estimator of the strength of the cross-correlation 
peak. During the course of this program — which spanned 
seven observing runs totaling 13 nights — we obtained 53 
distinct radial velocity measurements for 23 different pro- 
gram objects. All of these objects are faint, metal-poor 
members of our program clusters. Using the procedures 
described in Vogt et al. (1995) and Cote et al. (1999), we 
then derived empirical estimates for our radial velocity un- 
certainties. Specifically, we assume that the uncertainty of 
any radial velocity measurement, e{vr), can be expressed 



i{Vr)^a/{l+RTD) 



(2) 



where a is a constant to be determined. For our sample 
of repeat velocity measurements, which has 30 degrees of 



freedom, we calculate a = 9.0 



+2.4 



1.6 for X 



1, where the 



quoted uncertainties refer to 90% confidence limits. 

Table ^ summarizes the results of our photometric and 
spectroscopic observations of Palomar 13. From left to 
right, this table records the name of each program star,^ 
distance from the cluster center, V magnitude, {B — V) 
color, HIRES exposure time, heliocentric Julian date, 
Tonry & Davis Rtd value, heliocentric radial velocity, 
and the weighted mean velocity. The final column gives 
the cluster membership probability, P(/i), taken from the 
proper motion survey of Siegel et al. (2001). These mem- 
bership probabilities are discussed in detail below. 

3. RESULTS 

3.1. Identification of Members and Non-members 

A crucial first step in studying the internal dynamics 
of Palomar 13 is the isolation of a sample of bona fide 

^Identifications are from Ortolani et al. (1985) for those stars beginning with either a "ORS" or "V" prefix; four additional stars (#910, 
911, 915 and 931) which were not included in the Ortolani et al. (1985) study were selected from images taken with COSMIC on the Palomar 
5m telescope. 

^Ortolani et al. (1985) reported that the cluster luminosity function is similar to that of M3, but reduced by a factor of ~ 1/60. 

''Two stars have been omitted: V2 (a known RR Lyrae; Ciatti et al. 1965) and #915. Siegel et al. (2001) were unable to measure proper 
motions for either of these two objects, but both stars are almost certainly members based on their radial velocities, their location in CMDs, and 
their proximity to the center of Palomar 13 (see columns 2 of Table M and Figure hi). Nevertheless, V2 has been excluded from the calculation 
of the observed velocity dispersion since atmospheric pulsations may infiate the measured dispersion (e.g., the peak-to-peak radial velocity 
variations of RR Lyrae stars in M92 are Avr = 50 — 60 km s~^; Cohen 1992). 



members. At a Galactic latitude of 5 = — 42?7, Palomar 
13 is located well below the Galactic plane, but its very 
sparse giant branch® means that contamination from fore- 
ground disk and halo stars may be non-negligible. As a 
further complication, there is only limited radial velocity 
separation with the disk field star population due to Palo- 
mar 13's low systemic velocity (see §3.5). These difficulties 
are ameliorated by the high precision of our radial veloc- 
ity measurements, our new CMD, and the availability of 
proper motions from Siegel et al. (2001). 

In Figure |^, we plot the mean radial velocities of our 
program stars, (vr), against membership probability, P(/x), 
from Siegel et al. (2001).^ The sharp spike at (vr) ^ 25 
km s~^ identifies Palomar 13; the vertical dashed line indi- 
cates our best estimate for its systemic velocity (see §3.5). 
In general, there is good agreement between the member- 
ship classifications based on the radial velocities and those 
based on the proper motions. The three exceptions are 
ORS-23, ORS-32 and ORS-118. In the first case, the mea- 
sured radial velocity of Vr — 18.57 ±0.82 km s~^ for ORS- 
23 provides no compelling evidence either for, or against, 
membership. However, Siegel et al. (2001) find P{fi) = 
0% for this star, a conclusion which is supported by the 
location of this star in the CMD: i.e., the star has a very 
red color of {B — V) — 1.50 mag and is located ~ 0.8 mag 
blueward of Palomar 13's RGB. We conclude that ORS- 
23 is a likely field star. For ORS-32, Siegel et al. (2001) 
report an intermediate probability of P(/i) — 35%. Our 
mean radial velocity, Vr = 19.68 ± 0.26 km s~^, is based 
on four measurements at three different epochs. The resid- 
ual with respect to the systemic velocity of Palomar 13 is 
4.5 km s~^, and the star is located i? ~ 1' from the core 
of Palomar 13. We believe that these facts, coupled with 
the position of this star in the CMD (particularly in the 
V, V ^ I plane) , marginally favor the interpretation that 
it is a member of Palomar 13. Nevertheless, a spectro- 
scopic measurement of the metallicity of this star would 
be desirable. 

The situation for ORS-118 is somewhat ambiguous. Al- 
though Siegel et al. (2001) report P(/i) = 0% for this star, 
its location in the CMD strongly suggests that it is a true 
member, since it lies precisely on the RGB of the best-fit 
isochrone (see §3.3). The measured radial velocity of this 
star, Vr = 24.92 ± 0.21 km s~^, is indistinguishable from 
the systemic velocity of Palomar 13, and is based on three 
independent measurements which agree internally to bet- 
ter than 0.75 km s~^. We further note that an abundance 
analysis for ORS-118 (described below) yields [Fc/H] = 
— 1.98 ± 0.31 dex. This metallicity points to a halo ori- 
gin, and is consistent with the metallicity of Palomar 13 
deduced from isochrone fitting (see §3.3). 

To investigate the possibility that ORS-118 is simply 
an interloping foreground star which happens to have the 
same radial velocity as the Palomar 13, we approximate 
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the line-of-sight velocity distributions expected for disk 
and halo stars in this direction with the relation 

P{vr) oc exp [-{{vr) " (d,,o))V2ct'] (3) 

where 

(Vrfi) = -19.5 cos (4) 

is the mean radial velocity of disk stars along this line of 
sight as a result of the solar motion towards I = 56° and 
b = 23° (Mihalas & Binney 1981). The angle between this 
point on the celestial sphere and Palomar 13 is denoted by 
(j)d- Similarly, the mean radial velocity of halo stars along 
this line of sight is taken to be 

{Vr,o) = —220 COS (5) 

where (j)h is the angle subtended by Palomar 13 and {I, b) 
— (90°, 0°). We further approximate the local velocity 
dispersion of disk stars as 

+ (6) 

which, from the last row of Table 1 of Dehnen & Binney 
(1998), yields cr ~ 29 km s^^. For the halo, we assume an 
isotropic velocity ellipsoid (i.e., <Ju — o'v = o'z) and adopt 
a ~ 124 km s'^ from Cote (1999). 

The resulting probability distributions for disk and halo 
stars are illustrated by the dotted lines in Figure |[ Al- 
though the absolute normalization of the two curves is 
arbitrary, their relative normalization has been adjusted 
to match the expected number of disk and halo stars with 
16.5 < < 20 along this line of sight (determined from the 
IAS Galaxy model; BahcaU & Soneira 1980). If ORS-118 
is indeed a halo field star, then it lies 1.5a from the mean 
of the field distribution; the probability of such an occur- 
rence is 13%. Moreover, in the range 16.5 <V< 17.5, the 
IAS Galaxy model predicts only a single halo field star 
in our field, and just 8% of halo stars in this magnitude 
range are expected to have {B — V) > 0.93 (i.e., the color 
of ORS-118). We conclude that ORS-118 is likely to be a 
bona fide member of Palomar 13. The origin of the large 
proper motion residual measured by Siegel et al. (2001) 
is unclear, particularly since an inspection of our CCD 
images revealed no close companions. 

3.2. Metallicity 

Exposure times for our HIRES observations were cho- 
sen to yield the minimum signal-to-noise needed to derive 
radial velocities with a precision e{vr) <^ 1 km s~^. For 
ORS-118, the brightest star in our sample which has a ra- 
dial velocity consistent with membership in Palomar 13, 
we combined the two spectra obtained during the July 
1999 observing run. The signal-to-noise of ratio of this co- 
added spectrum — S/N ~ 25 per resolution clement — is 
just adequate for abundance analysis. 

Our analysis used 28 Fe I lines with oscillator strength 
values adopted from McWilliam et al. (1995), Kurucz 
model atmospheres and solar abundances from Anders & 
Grevesse (1989). Computations were made with the most 
recent version of the LTE line analysis code MOOG (Sne- 
den 1973). The effective temperature, Te, was first de- 
termined using an initial estimate from the photometry 



and then checked through excitation equilibrium of the 
Fe I lines whose equivalent widths were measured using 
standard IRAF routines. The surface gravity, log 17, was 
derived from the ionization equilibrium of Fe I and Fe II 
lines (using four Fe II lines in the available spectral re- 
gion). Once the effective temperature and surface gravity 
were determined, abundances for the Fe I lines were cal- 
culated by iteratively varying the microturbulent velocity 
until the best fit in the diagram of Fe I abundance versus 
equivalent width was achieved. Our best-fit parameters 
for ORS-118 are = 4700 K, logg = 1.7, and [Fe/H] 
= —1.98 ± 0.31 dex, where the rather large uncertainty 
on the measured metallicity is a reflection of the modest 
signal-to-noise of our co-added spectrum. 

Previous estimates for the metallicity of Palomar 13 in- 
clude [Fe/H] = -1.9 ± 0.4 (Canterna & Schommer 1978), 
-1.67 ±0.15 (Zinn & Diaz 1982) and -1.9 ±0.1 dex (Friel 
et al. 1982). The measurement of Canterna & Schommer 
(1978) comes from Washington photometry of individual 
stars, while the latter two measurements are based on low- 
resolution spectroscopy. Our determination of the cluster 
metallicity is slightly lower than, but still consistent with, 
these previous estimates. In the following section, we show 
that this metallicity is also consistent with that found from 
isochrone fitting of the CMD. 

3.3. Color Magnitude Diagram 

Figure ^ shows BV and VI CMDs for Palomar 13 con- 
structed from our LRIS photometry. As noted by previous 
investigators, the RGB is very sparsely populated. This 
figure (and the left panel of Figure ||) also confirms the 
finding of Siegel et al. (2001) that Palomar 13 contains 
a population of blue stragglers; likwise, there may be ev- 
idence for a "second sequence" which runs parallel to the 
upper main sequence. Both this latter sequence (which 
is traditionally interpreted as the signature of unresolved 
binary stars; e.g., Romani & Weinberg 1991), and the pres- 
ence of blue stragglers (a fraction of which are likely to be 
W UMa systems; Mateo et al. 1990) provide some evi- 
dence for a population of binary stars. If this is indeed 
the case, then it would have important implications for 
modeling the observed velocity dispersion profile. In §5.1, 
we shall return to the question of binary stars and their 
possible effect on the observed velocity dispersion. 

From the DIRBE maps of Schlegel, Finkbcincr & Davis 
(1998), we find the reddening in the direction of Palomar 

13 to be E{B - V) = 0.11 mag. The sohd curve in the 
left panel of Figure ^ shows a 14 Gyr isochrone from Berg- 
busch & VandenBerg (1992) with [Fe/H] = -1.78 dex, 
shifted by (m — M)v = 17.27 mag. We stress that the 
sparsely populated red giant and horizontal branches pre- 
clude a precise age determination, and simply note that a 

14 Gyr isochrone provides a reasonable match to the ob- 
served fiducial sequences. Assuming Ay = 3.1E{B — V), 
we find a true distance modulus of (m — M)o = 16.93±0.10 
mag, which corresponds to a distance of D — 24.31}'^ 
kpc. This value is in good agreement with the distance 
of D = 24.8 estimated by Siegel et al. (2001) from the 
average magnitude of a small sample of cluster RR Lyrae 
and horizontal branch stars. Based on the combination of 
isochrone fitting to the cluster CMD and the spectroscopic 
analysis of ORS-118, we adopt [Fe/H]=-1.9 ± 0.2 dex as 
our best estimate for the metallicity of Palomar 13. 
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Ortolani et al. (1985) found an absolute magnitude of 
My ~ —3.4 mag for Palomar 13 by integrating the lu- 
minosity function. For a distance oi D = 24.3, we find 
an integrated luminosity of 2.4 x 10'^ Ly Q {My — —3.6) 
within the tidal radius suggested by Siegel et al. (2001), 
rt ~ 188" (although, as discussed in §3.4, we see no evi- 
dence for a tidal cutoff at this, or any other, radius). This 
estimate excludes stars with [B — V) < 1.3 mag and con- 
firmed non-members. If the integration is extended over 
the entire LRIS field of view (as would be appropriate 
for the larger tidal radius found in §4), the luminosity in- 
creases to 3.1 X 10'^ iy,0 {My = —3.9). Our corresponding 
best estimate for the absolute magnitude of Palomar 13 is 
My = —3.8 mag (see also §4-5). In the most recent version 
of the Harris (1996) catalog of Galactic globular clusters, 
only three objects — E 3, Palomar 1 and AM 4 — have 
lower luminosities.® The integrated, de-reddened color is 
{B — V)q — 0.65 ± 0.04 mag, roughly consistent with that 
expected for a globular cluster of metallicity [Fe/H] ~ —1.9 
dex (Durrell et al. 1996). 

3.4. Surface Density and Surface Brightness Profiles 

While there are surprisingly few published studies of 
the structural properties of Palomar 13, previous results 
have consistently pointed to a compact, faint cluster of 
low concentration. For instance, Webbink (1985) quotes 
a core radius of Tc — 23" and a concentration index of 
c = log (rt/rc) = 0.9, while Trager et al. (1995) report 
Tc = 29" and c — 0.66. However, both of these studies 
relied primarily on the star counts of Kinman & Rosino 
(1962) which were obtained by viewing the projection of 
a photographic plate through a mask fixed on the screen 
of a Sartorious astrophotometer. Recently, Siegel et al. 
(2001) adopted c = 0.7 and measured Tc = 39" for Palo- 
mar 13 using a sample of 119 members selected on the 
basis of their proper motions, but noted that this model 
provides a poor description of the measured profile at large 
radii: i.e., 14 of their proper motion members (12% of the 
sample) lay beyond the canonical tidal radius. This fact, 
coupled with Palomar 13's low luminosity and eccentric 
orbit, led Siegel et al. (2001) to suggest that this object is 
in the final stages of destruction. 

We have used our Keck and CFHT imaging to mea- 
sure improved structural parameters for Palomar 13. The 
point-source catalog from our Keck images consists of 841 
objects distributed over an area of 42.8 arcmin^. Since 
the Palomar 13 field contains a significant number of obvi- 
ous foreground/background objects — see Figure |— the 
Keck point-source catalog was trimmed to exclude objects 
with 0.0 < {B — V) < 1.3 mag. Similarly, sources which 
did not fall in the range 0.15 < {V — r) < 0.7 were culled 
from the CFHT catalog. The point-source catalogs were 
further restricted to objects brighter than V — 24: (Keck) 
and 23.5 (CFHT), in order to ensure photometric com- 
pleteness at all radii. This leaves a total of 629 and 551 
objects in the Keck and CFHT databases, respectively. 

Figure |^ shows the surface density profiles determined 
from these catalogs. Observations from Keck are shown 
as open squares, while open circles denote CFHT data 
points. The latter have been scaled upward by a factor 
(629/551) ~ 1.14 to account for the difference in sample 

^According to Harris (1996), the absolute magnitude of Terzan 1 is 
have recently reported an upward revision to My ~ —5.4 for this clusi 



size. There is good agreement between these two indepen- 
dent profiles. Given the larger radial extent of the CFHT 
data, they provide stronger constraints on the asymptotic 
behavior of the surface density profile — a particularly 
important issue in the present case since the Keck surface 
density profile continues to decline out to the edge of the 
LRIS field of view. The arrow in Figure || shows our best 
estimate for the background surface density, E;, = 0.68 
stars arcmin"^, based on the scaled CFHT observations. 
For comparison, the surface density of Galactic stars in the 
appropriate color and magnitude intervals in the direction 
of Palomar 13 predicted by Institute for Advanced Study 
Galaxy model (Bahcall & Soneira 1980; Bahcall 1986) fall 
in the range 0.64 ^ Sf, ^ 0.95 stars arcmin^^, consis- 
tent with our empirical estimate. The filled symbols in 
this figure show the background-subtracted surface den- 
sity profile. 

For comparison, the open triangles in Figure |^ shown 
the surface density profile of Siegel et al. (2001), taken di- 
rectly from a digitized image of their Figure 7. The Siegel 
et al. data points have been scaled upwards by the ratio 
of the numbers of stars in the two studies: (629/119) ~ 
5.29. In the inner regions of Palomar 13, the agreement 
between the two profiles is excellent. Beyond a distance 
of i? = 2 — 3', though, the profiles differ, with our sur- 
face densities being higher by a factor of ^ five (although, 
to within the rather large uncertainties, the slopes of the 
two profiles are in agreement). The origin of this differ- 
ence is unclear. The profile of Siegel et al. (2001) has the 
obvious advantage of being proper motion selected: i.e., 
objects having membership probabilities of P(/i) < 50% 
have been rejected on a star-by-star basis. In principal, 
the use of proper motions to cull interlopers should be su- 
perior to the statistical approach used here to estimate the 
background level. On the other hand, Siegel et al. (2001) 
note that most of their plates extend to only V ^ 21 (i.e., 
the approximate location of the main-sequence turnoff) 
whereas each of our two independent sets of CCD pho- 
tometry reach several magnitudes below this level, and 
with small photometric errors. If Palomar 13 is experi- 
encing significant mass loss through evaporation and tidal 
stripping (see §5.2), then it is the lowest mass stars which 
are expected to have largest evaporation rates (i.e., they 
should preferentially inhibit the cluster envelope as a result 
of mass segregation). In such a case, our deeper star counts 
compared to Siegel et al. (2001) might lead to a more pro- 
nounced surface density excess at large radii. In any event, 
we caution that the five outermost bins in the Siegel et al. 
(2001) profile (spanning the range 2' <^ R <^ 10') contain 
just 15 stars, so it is unclear of the difference between the 
two studies is truly significant. 

With these caveats in mind, we point out a striking fea- 
ture of Figure I — the near power-law behavior of the 
surface density profile. Indeed, a single power-law, E(i?) 
oc R ' , with 7 — 1.8 ± 0.2 provides an adequate descrip- 
tion of the surface density profile for R ^ 0'.3. The dotted 
curve in Figure ^ shows the King-Michie model proposed 
by Siegel et al. (2001). This model, which has a somewhat 
arbitrary concentration parameter of c = 0.7 and a core 
radius of Vc = 39", provides adequate description of the ob- 

^ = —3.3; i.e., fainter than Palomar 13. However, Idiart et al. (2002) 



COTE ET AL. 



7 



served surface density profile over the range 0.'3 <i -R ^ 1', 
but underestimates the central surface density. More im- 
portantly, if this model accurately reflects the true density 
profile of Palomar 13, there is an even more dramatic "ex- 
cess" of stars above the model predictions beyond R 1'. 
While it is certainly true that the exact shape of the sur- 
face density profile depends on the adopted background 
level (particularly at large radii) this excess is unlikely to 
be an artifact of the background subtraction. For instance, 
achieving an acceptable fit to a King model with c <^ 1 for 
R <^2' would require the background level to be increased 
by roughly an order of magnitude. 

The dashed curve in Figure || shows the best-fit King- 
Michie model, which will be discussed in detail in §4. As 
shown there, the King-Michie which best fits the observed 
surface density and surface brightness profiles of Palomar 
13 has a tidal radius of rt — 26±6' (182±41 pc). Since this 
is more than twice the radial extent of our CFHT imaging, 
the adopted background level should be viewed with some 
caution. However, we note that estimating the background 
surface density from star counts made within rt should re- 
sult in an overestimate of the true background level and, 
hence, to an underestimate of rt- We are therefore left with 
two possible interpretations of the surface density profile of 
Palomar 13: either this object contains a significant popu- 
lation of "extra-tidal" stars which have been caught in the 
act of evaporation (as suggested by Siegel et al. 2001), 
or its concentration and spatial extent have been greatly 
underestimated in previous studies. 

We now turn our attention to the properties of the stars 
associated with this surface density excess. This feature 
is evident in both the CFHT and LRIS surface density 
profiles, but is there any evidence that this population is 
truly comprised of Palomar 13 members? Figure ^ shows 
instrumental CMDs for the Palomar 13 field derived from 
our CFHT photometry, along with the color and magni- 
tude limits used to cull probable non-members from the 
sample (i.e., the dotted regions). The panel on the left 
shows the CMD for all stars within the Siegel et al. (2001) 
tidal radius. The middle panel shows the CMD for all ob- 
jects beyond this radius, while the right panel shows those 
stars located at more than twice this distance. The dashed 
curves in the final two panels show the Palomar 13 fiducial 
sequence. A significant number of probable members with 
R > rt are evident, providing additional evidence that the 
stars belonging to the excess in the surface density profile 
are indeed associated with Palomar 13. 

Figure shows the surface brightness profile for Palomar 
13 measured from our Keck and CFHT images. Profiles 
having been calculated independently using both the Keck 
and CFHT photometry. Surface photometry is shown as 
the pentagons (Keck) and triangles (CFHT), while the re- 
sults from star counts are indicated by the squares (Keck) 
and circles (CFHT). A direct measurement of the back- 
ground surface brightness for the Keck and CFHT im- 
ages yields = 21.3 and 21.5 mag arcsec"^, respectively 
— consistent with expectations for dark sky conditions 
at Mauna Kea during 1999 (Krisciunas 1997). In deriv- 
ing a surface brightness profile from the star counts dis- 
cussed above, we have simply multiplied the surface densi- 
ties shown in Figure |^ by a constant factor, rj — Ly/ritot, 
where Ly is the integrated luminosity of Palomar 13 and 



ntot are the number of stars used to derive each of the 
Keck and CFHT surface density profiles. This approach 
suppresses large fluctuations in the surface brightness pro- 
file caused by the small number of red giants in this low 
luminosity object. 

How does our surface brightness profile compare with 
previous findings? Unfortunately, Ortolani et al. (1985) 
and Siegel et al. (2001) each presented surface density pro- 
files for Palomar 13, but neither reported surface bright- 
ness measurements. Webbink (1985) quotes a central sur- 
face brightness of nv = 22.41 mag arcsec"^ based on the 
star counts of Ortolani et al. (1985) and the aperture 
photometry of Racine (1975). Trager et al. (1995) found 
a similar core radius and concentration (i.e., Tc = 29" and 
c — 0.66) but a much different central surface brightness: 
/iy — 24.31 mag arcsec"^. The overall shape of the Trager 
et al. (1995) surface brightness profile, shown by the open 
stars in Figure |7[ is in reasonable agreement with the pro- 
file found here, but differs systematically by a factor of ~ 
3.4 in surface brightness (in the sense that their measure- 
ments are lower). Their profile consists primarily of star 
counts Kinman & Rosino (1962), normalized to the aper- 
ture photometry of Racine (1975). We believe that the 
normalization given here is to be preferred since two inde- 
pendent measurements (from Keck and CFHT) are in close 
agreement. Moreover, an integration of the best-fit King- 
Michie model of Trager et al. (1995) yields an absolute 
magnitude for Palomar 13 of My — —1.25 mag. This is ^ 
ten times fainter than the value of My — —3.6 reported by 
Ortolani et al. (1985) and the value of My — —3.8 found 
here, both of which were obtained by direct integration of 
the luminosity function. 

The various curves shown in Figure |^ represent different 
models for the light distribution in Palomar 13. We will 
discuss these models in detail below (in §4-5), but before 
doing so, we shall use the radial velocities presented in Ta- 
ble 2 to investigate the internal kinematics of Palomar 13. 

3.5. Internal Kinematics 

Using our sample of 21 members (V2 excluded), we cal- 
culate a systemic velocity of (wr}s — 24.1 ± 0.5 km s~^ 
using the maximum likelihood estimator of Pryor & Mey- 
lan (1993). The corresponding estimate for the the pro- 
jected intrinsic velocity dispersion using this same sample 
is CTp = 2.2 ± 0.4 km s'^. Omitting ORS-118 from the 
sample (see §3.1) has essentially no effect on the derived 
parameters: {vr) = 24.0 ± 0.5 km s~^ and Up = 2.3 ± 0.4 
km s~^. Likewise, if ORS-32 is omitted from the sam- 
ple, the best-fit parameters do not change significantly: 
(vr) = 24.3 ± 0.5 km s~i and CTp = 2.0 ± 0.4 km s'^. 

In the upper panel of Figure |[ we show the radial ve- 
locities for these stars plotted against distance from the 
center of Palomar 13. The horizontal dashed line shows 
the systemic velocity, {vr)s = 24.1 km s~^, while the ver- 
tical arrows aX R = V denote the radius where the possible 
excess in the surface density profile begins. It is perhaps 
noteworthy that the stars with the largest velocity resid- 
uals are located at, or beyond, this radius. We note that 
there is a tendency for stars with the most discrepant ve- 
locities to also have the smallest uncertainties — this is a 
consequence of the fact that these stars were observed on 
multiple observing runs to confirm the first-epoch veloci- 
ties and to search for any evidence of binarity. As with the 



8 



PALOMAR 13 



full sample, none of these stars shows compelling evidence 
for radial velocity variations. 

In the lower panel of Figure |[ we plot the velocity dis- 
persion profile for Palomar 13. The open circle denotes the 
dispersion found using our complete sample of 21 mem- 
bers, plotted at their mean radius: {R) — 0.'84 (5.9 pc). 
The filled circles show the results of dividing the sample 
into two radial bins containing 11 and 10 stars, respec- 
tively. In this case, the dispersions are 1.0±0.4 km s~^ at 
(R) = 0.'36 (2.6 pc) and 2.6±0.6 km s"! at (i?) = lf36 
(9.6 pc). Although there may be some evidence for a ris- 
ing velocity dispersion profile, we caution that this result 
is based on a small sample, and that the uncertainties in 
the measured dispersions are appreciable. 

How does the measured velocity dispersion for Palomar 
13 compare with expectations? Globular clusters are well 
described by King-Michie models, in which the central ve- 
locity dispersion (in km s^^) may be approximated by 

CTp - 0.003TyrcMy,o (7) 

where Ty is the V^-band mass-to-light ratio, Vc is the core 
radius in pc, and ^v,o is the central surface brightness is 
units of Lv,Q pc~^ (e.g., Richstone & Tremaine 1986; Ma- 
teo et al. 1991). Based on the King-Michie model fits 
to the observed surface brightness profile presented in §3.4 
and Table H, the expected velocity dispersion of Palomar 13 
is (Tp.o — 0.6 km s^^ for an assumed mass-to-light ratio of 
Ty = 2. Thus, our measured velocity dispersion exceeds 
the expected value by roughly a factor of four. The im- 
plications of this finding are presented below, along with 
a more detailed analysis of the internal kinematics Palo- 
mar 13. 

4. MASS-TO-LIGHT RATIO 

As discussed in §3.4, our new surface density profile for 
Palomar 13 suggests that either a significant fraction of 
its stars are located beyond the ostensible tidal radius, or 
the object is more spatially extended than previously be- 
lieved. We now examine its internal kinematics within the 
context of these two scenarios. 

The short-dashed curve in Figure ^ shows the King- 
Michie which best fits the observed surface brightness pro- 
file for Palomar 13. The parameters for this model are = 
14±2", c = 2.0±0.15 and fiv,o = 22.541°;?? mag arcsec-^. 
(These and other model parameters are summarized in Ta- 
ble ^, along with various observed and derived properties 
of Palomar 13.) According to this model, Palomar 13 has 
both a higher concentration and a greater spatial extent 
(rt ~ 180 pc) than suggested by previous studies. Indeed, 
this concentration index is considerably larger than that 
expected for a globular cluster of this luminosity or surface 
brightness (see, e.g., Djorgovski 1991). Integration of this 
profile yields a total luminosity of Ly — (2.8 ± 0.4) x lO'^ 
Lv,Q- The corresponding velocity dispersion profile for 
an assumed mass-to-light ratio of Ty = 2 is shown by 
the lower dashed curve in the bottom panel of Figure @ 
(i.e., curve 1-A). In agreement with the naive analysis 
presented above, we find the measured velocity dispersion 
to exceed that expected for a normal globular cluster of 
this concentration and luminosity. A maximum likelihood 
solution for the best-fit scale velocity of this same King- 
Michie model gives the upper dashed curve shown in this 



figure (i.e., curve 1-B). This model has a central veloc- 
ity dispersion of o — 2.5 km s^^, and a total mass of 
M = (1.10ln;36) Mq. The mass-to-light ratio is then 
Ty = 40l^y (1(7 uncertainties). 

Alternatively, we may use the King-Michie model sug- 
gested by Siegel et al. (2001) to estimate the mass and 
mass-to-light ratio of Palomar 13. This model, which is 
shown by the dotted curve in Figure |^, yields a total lumi- 
nosity of Ly ~ 1.2 X 10"^ Lv,Q- A maximum likelihood scal- 
ing of the model velocity dispersion profile to the measured 
velocities produces the dotted curve in the lower panel of 
Figure^ (i.e., curve 2). For this model, which has a central 
velocity dispersion of cr^.o — 2.8 km s^^, the mass within 
the tidal radius {n = 188" or 23 pc) is M = x 10"^ 

Mq. The corresponding mass-to-light ratio is Ty = 48li[4 
where the quoted Icr uncertainties refer only to the errors 
on the derived mass. Thus, irrespective of the adopted sur- 
face brightness profile, it is not possible to reproduce the 
observed velocity dispersion with a mass-to-light ratio that 
is typical of globular clusters (e.g., (Ty) — 2.3±1.1; Pryor 
& Meylan 1993). 

5. IMPLICATIONS 

What are the implications of this finding? We begin 
by examining — and tentatively rejecting — the possibil- 
ity that the observed velocity dispersion has been inflated 
by a "velocity jitter" in the atmospheres of the program 
stars, or by a population of spectroscopic binary stars. We 
then examine other scenarios which may explain the ob- 
served properties of Palomar 13 including tidal disruption, 
a mass function which is heavily skewed to low-mass stars 
or heavy stellar remnants, modified Newtonian dynamics, 
and/or the presence of a massive dark halo. It is worth 
bearing in mind that these various possibilities may not be 
mutually exclusive (i.e., a very high binary fraction com- 
bined with tidal heating and mass loss). 

5.1. Velocity Jitter and/ or Binary Stars? 

The value of dp — 2.2 ± 0.4 km s^^ reported in §3.5 
refers to the intrinsic velocity dispersion for Palomar 13: 
i.e., the contribution from measurement uncertainties has 
been removed (see, Pryor & Meylan 1993). Our veloc- 
ity uncertainties have been calculated empirically from 53 
distinct radial velocities of faint metal-poor stars (accumu- 
lated during multiple observing runs) so we are confident 
that the measured dispersion is not the result of underes- 
timated velocity uncertainties. Such an explanation would 
require a ~ 20 km s~^ (see Equation H), whereas we mea- 
sure a = 9.0 km s^^ and find a < 12.8 km s'^ at 99% 
confidence. 

Bright red giants in globular clusters often show evi- 
dence for radial velocity variations in the range 2-5 km 
s~i (Gunn & GrifRn 1979; Mayor et al. 1984). This so- 
called velocity "jitter" is thought to arise from pulsational 
motions in the atmospheres of these evolved stars. It is 
unlikely that our measured velocity dispersion is the re- 
sult of such velocity variations since long-term monitoring 
of red giants in nearby globular clusters has established 
that the magnitude of this jitter shows a strong luminos- 
ity dependence: i.e., it is significant at the 1 km s~^ level 
only for stars which are located within ~ 1 mag of the 
RGB tip. Our radial velocity sample includes no stars 
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this bright. In fact, our most luminous member (ORS- 
118) has My ~ —0.3 mag — a full 2.3 mag below the tip 
of the RGB. 

Might our estimate for the velocity dispersion be the 
inflated by "dynamically-hard" spectroscopic binary stars 
residing in the cluster? There are reasons to believe that 
the binary star fraction in Palomar 13 may be appreciable. 
First, the CMD of Palomar 13 shows both a sizeable blue 
straggler population and a possible "second sequence" — 
features which are usually interpreted as signatures of un- 
resolved binaries (Mateo et al. 1990; Romani & Weinberg 
1991). Second, on purely theoretical grounds, the low stel- 
lar density and velocity dispersion of Palomar 13 (for an 
assumed "normal" mass-to-light ratio) suggest that even 
wide binaries might have survived disruption over a Hub- 
ble time. For instance, for a central velocity dispersion of 
(Tp^o = 1 km s~^, a pair of equal-mass 0.8Mq stars would 
be disrupted by encounters with other members only if 
separated by more than Ocrit ~ 1,400 AU (Hills 1984). 
Such binaries, if on circular orbits, would have periods of 
P ~ 4 X lO'' yr and orbital velocities of ~ 1 km s~^: i.e., 
modest, but potentially important, given the low disper- 
sion expected for Palomar 13. 

To explore the possible effects of binary stars on our 
measured velocity dispersion, we have carried out Monte- 
Carlo simulations similar to those described in Pryor, 
Latham & Hazen (1988) and Cote et al. (1996). In brief, 
we have generated simulated radial velocity measurements 
having the identical precision and temporal spacing as the 
actual dataset. For each simulation, we take the intrin- 
sic velocity dispersion to that predicted by a King-Model 
having Ty ~ 2, and calculate the contribution to the 
measured velocity dispersion of a randomly drawn per- 
centage, /b, of unresolved spectroscopic binary stars. Bi- 
naries are assumed to have circular orbits. The orbital 
periods in years, P, and mass ratios, q = M2/M1, are 
drawn at random over the ranges, — 1.0 < logP < 4.0 and 
—0.6 < logg < 0.0. For each /;,, we generate 1000 arti- 
ficial datasets and calculate the mean cluster dispersion 
in precisely the same manner as the actual for the actual 
data; the percentage of binaries is varied until the median 
dispersion found in 1000 simulated datasets is equal to the 
value obtained using the actual observations. These simu- 
lation suggest that a value of fb — 30% produces a median 
velocity dispersion which best matches the observed value 
of cTp = 2.2 km s^^. If it is assumed that the binaries 
do not have circular orbits, but instead have a "thermal" 
distribution of eccentricites, then the percentage of bina- 
ries needed to explain the measured dispersion would be 
roughly twice as large (see Pryor et al. 1988; Cote et al. 
1996). 

By definition, the velocity dispersion profile predicted 
by these simulations differs from that of the input King- 
Michie model only in normalization since we have made 
the first-order assumption that the binaries have the same 
density distribution as the single stars. However, if the 
binaries are instead assumed to be more centrally con- 
centrated than the other stars — perhaps as a result of 
mass segregation — then the expected velocity dispersion 
profile would show a steeper decline. Needless to say, a 
rising velocity dispersion profile would prove difficult to 
reconcile with this scenario. We find weak evidence for 
an outward increase in the observed velocity dispersion 



of Palomar 13 but given the small number of measured 
velocities, this feature clearly needs to be confirmed by 
additional observations. As Figure || demonstrates, there 
are many promising candidates beyond R ^ 2' (i.e., the 
radii of the most distant members in the current radial 
velocity sample). It would be of considerable interest to 
measure velocities for these stars in order to determine if 
they are members of Palomar 13 and, if so, to see if the 
velocity dispersion profile rises with increasing radius. 

Likewise, it would be desirable to search directly for ve- 
locity variations associated with possible orbital motions 
of these stars. Unfortunately, the long periods of the pu- 
tative binaries would render any such search hopeless; a 
more promising approach would be to measure fb by mod- 
eling the "second sequence" measured from high-resolution 
T/ST imaging (e.g., Rubenstein & Bailyn 1997). 

Finally, we note that spectroscopic binaries may offer a 
viable explanation for the observed high velocity disper- 
sion, but they cannot on their own explain the anomalous 
surface density profile. 

5.2. Disruption? 

An underlying assumption of the dynamical analysis 
presented in §4 is that of virial equilibrium. In light of 
the large velocity dispersion measured for Palomar 13, it 
is prudent to investigate the validity of this assumption, 
which has recently been called into question for many of 
the dwarf satellites and globular clusters in the Galactic 
halo (Majewski et al. 2001). Simply put, could Palomar 
13 be in the process of disruption? This possibility was 
advocated by Siegel et al. (2001), who combined their 
proper motion measurements with the systemic radial ve- 
locity presented here to derive an orbit for Palomar 13. 
They concluded that Palomar 13 is just ^ 70 Myr past 
perigalacticon {Rp ~ 11.2 kpc) and suggested that it has 
undergone catastrophic tidal heating during its recent ap- 
proach. Since N-body simulations of tidal tails in globular 
clusters (Combes, Leon & Meylan 1999) show that the 
time taken for evaporating stars to diffuse along the clus- 
ter orbit can exceed this timescale, it is worth examining 
this possibility further. 

Siegel et al. (2001) noted the presence of a population 
of "extra-tidal" stars in their surface density profile. We 
too find a possible excess in the Palomar 13 surface den- 
sity profile, beginning at i? ^ 1'. This feature is evident in 
both the CFHT and LRIS surface density profiles, and the 
CFHT photometry suggests that many of these stars have 
magnitudes and colors which are consistent with those ex- 
pected for Palomar 13 members (see §3.4). Thus, we con- 
firm the existence of this feature; the issue is whether the 
stars associated with it are truly evaporating from Palo- 
mar 13, or if they are instead bound members. 

As shown in Figure ||, the surface density profile 
for Palomar 13 shows a roughly power-law decline for 
R ^ 0f3. This behavior is most pronounced beyond 
R ^1', where there may be some evidence for a "break" 
in the profile. Johnston et al. (1999; 2001) have car- 
ried out N-body simulations of satellites orbiting in the 
Galactic tidal field; they find that the radial distribution 
of unbound stars in the outer regions of disrupting satel- 
lites should define a power-law profile — broadly similar 
to what we observe for Palomar 13. The dashed lines in 
Figure |^ show power-law profiles, S(i?) oc R having 
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representative exponents of 7 = 0.5, 1.5 and 2.5. Over 
the range 1' <^ R <^ 4', we measure 7 = 1.8 ± 0.2 for Pafo- 
mar 13. This value is steeper than the formal value of 
7 ~ 1 estimated by Johnston et al. (1999) but is probably 
consistent with expectations given uncertainties in the nu- 
merical treatment (e.g., the initial density profile, viewing 
geometry and orbital phase; Johnston et al. 2001). 

If this interpretation is correct, then a crude estimate 
for the mass loss rate may be made as follows. Johnston 
et al. (1999) present two simple expressions for the frac- 
tional mass loss rate using extra-tidal stars, 

{df/dt)i = cos{e)TrRtN,t/[iRcct - Rb)Ni,P] 
{df/dt)2 = cos{e)^{Rb)2^Rl/[NbP], 

where 9 is the angle subtended by our line of sight and 
the plane perpendicular to the satellite's motion, Rjj in 
the "break" radius where the surface density excess be- 
gins, Rxt is the radius where the profile of the extra-tidal 
component ceases to be well defined, Nj;t is the number of 
stars in the range Rb < R < Rxt, and P is the satellite's 
azimuthal period (taken to be 2.0 Gyr following Dinescu 
et al. 2001). Using the space velocity of Siegel et al. (2001) 
and the surface density profile shown in Figure o, we find 
{df/dt)i ~ 0.2 Gyr-i and {df/dt)2 0.15 Gyr^. There 
are significant uncertainties involved in these calculations 
but it is clear that, if this interpretation of the surface den- 
sity profile is correct, then both estimators point to enor- 
mous recent mass loss. The inferred mass loss rate is far 
higher than those estimated by Leon, Meylan & Combes 
(2000) from star counts in a sample of 20 Galactic globu- 
lar clusters — e.g., in the case of oj Centauri, Leon et al. 
(2000) find that only « 1% of the cluster's mass has been 
lost during the last tidal encounter. Taken at face value, 
the mass loss rates for Palomar 13 suggest that it has suf- 
fered severe heating, and will disrupt entirely within the 
next few Gyrs. 

On the other hand, there may be some difficulties with 
this interpretation. N-body simulations of the tidal heat- 
ing of satellites indicate that the tidally stripped mate- 
rial tends to diffuse along the orbits, forming leading and 
trailing streams. We have examined our star counts for 
evidence of such streams by measuring the mean surface 
density of stars in four quadrants over the radial range 
1' < i? < 2.'88. The inner limit corresponds to the the 
break radius, while the outer limit is the most distant 
radius at which our geometric coverage is still 100% com- 
plete. The results of this exercise are shown in Figure ^ 
We find only weak evidence that these stars are oriented 
along a preferred direction. The density maxima of the 
best-fit sinusoid shown in this figure are found at position 
angles of Op = 15 ± 11° and 195 ± 11°. The agreement 
between this orientation and those of the absolute proper 
motion vector, 8p — 48°, and the Galactic center direc- 
tion, 8p — 240°, is unexceptional: i.e., the differences in 
position angles are AOp = 33 ± 11° and 45 ± 11°, respec- 
tively. There is poor agreement with the vector perpen- 
dicular to the Galactic disk, which has 8p = 327°. If the 
peculiar surface density profile of Palomar 13 is indeed 
the result of tidal heating, then this lack of a preferred 
orientation for the putative extra-tidal stars is somewhat 
puzzling. It would be worthwhile investigating this issue 
using star counts which extend over a wider field, prefer- 



ably well beyond the modest range radial range (i? <^ 3') 
examined above. 

Piatek & Pryor (1995) have used N-body simulations 
to show that it is possible to heat the outer regions of 
satellites during close passages with the Galaxy. Such en- 
counters often have the effect of inducing ordered motions 
at large radii while leaving the central velocity dispersion 
unchanged. According to the simulations, the gradient 
in velocity is usually maximized along the vector which 
points to the Galactic center. This vector is shown by 
the long arrow in Figure |^; for comparison, the short ar- 
row shows the direction of the cluster motion as measured 
by Siegel et al. (2001), while the intermediate arrow de- 
notes the direction perpendicular to the Galactic plane. 
As shown in Figure |l^, we find no statistically significant 
evidence for "rotation" among our program stars, as might 
be expected for stars in the leading or trailing streams of 
a disrupting cluster (see also Figure |^ where the velocity 
residuals are illustrated on a star- by-star basis). 

Finally, we note that there is an additional (albeit in- 
direct) constraint on the mass and dynamical status of 
Palomar 13. The tidal radius of a satellite object orbiting 
in the Galactic potential is given by 

rt = i?p[M/MG(3 + e)]i/3 (9) 

where Rp is the distance from the Galactic center at peri- 
galacticon, Mq is the interior Galactic mass at this point, e 
is the orbital eccentricity of the satellite, and M is its mass 
(King 1962). If we assume a canonical mass-to-light for 
Palomar 13 (e.g., Ty = 2), we may use equation (9) to ob- 
tain a crude estimate for the expected tidal radius. Adopt- 
ing e = 0.76 and Rp — 11.2 kpc for Palomar 13 (Siegel 
et al. 2001), and assuming using AIoiRp) = Oc^Rp/G 
(Fich & Tremaine 1991) and 9^ = 220 km s"\ we find 
rt ~ 26 pc for a total luminosity of Ly = 3 x lO'^ Lv,q- 
This is remarkably close to the tidal radius of rt = 23 pc 
for the King-Michie model of Siegel et al. (2001), which is 
shown as the dotted curves in Figures |andg. This IS con- 
sistent with the view that the excess in the surface density 
profile is a population of extra-tidal stars. The alternative 
interpretation — that the much larger tidal radius implied 
by the best-fit King-Michie model in Figure ^ is evidence 
for dark matter in Palomar 13 — will be examined in §5.4 

5.3. Baryonic Dark Matter: Low-Mass Stars and/or 
Heavy Remnants? 

Taken at face value, the measured velocity dispersion 
for Palomar 13 implies a total mass of M ~ (LIIq ^) x 
10^ Mq, far greater than that expected on the basis of 
the observed luminosity: i.e., M = TyLy ~ 6 x 10'^ Mq. 
The obvious baryonic candidates for this additional mass 
would be low-mass main-sequence stars, or heavy stellar 
remnants such as white dwarfs and neutron stars. Both va- 
rieties of heavy remnants are known to populate globular 
clusters (see Hut et al. 1992) and the contribution of white 
dwarfs to the overall mass budget is thought to increase 
as clusters evaporate (e.g., Vesperini & Heggie 1997). In 
this case, the low-luminosity of Palomar 13 might be in- 
terpreted as indirect evidence for an advanced stage of 
dynamical evolution. However, if one accepts the view 
that Palomar 13 is in dynamical equilibrium (see above). 
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then this would require ~ 95% of the total mass to re- 
side in heavy remnants. This percentage is far higher 
than that found in nearby, well-studied clusters: e.g., Mey- 
lan (1987; 1989) finds the percentage of heavy remnants 
(mainly white dwarfs) in 47 Tucanae and lu Centauri to 
be 20 ± 5% and 25 ± 15%, respectively (though the dy- 
namical mass found here for Palomar 13 is only half that 
contained in heavy remnants in 47 Tuc, by virue of the 
latter 's much larger mass). Unfortunately, a direct search 
for white dwarfs in Palomar 13 would be extremely chal- 
lenging since even the brightest such objects are expected 
to have ^ 27 mag at this distance (Richer et al. 1997)^ 
On the other hand, a determination of the main- 
sequence luminosity function down to this level (corre- 
sponding to My = 10 mag and Ty — 25; Silvestri et al. 
1998) is well within the reach of HST, and could put in- 
teresting constraints on the mass contained in low-mass, 
main-sequence stars. However, as shown in Figure ^ the 
available data suggest that the main-sequence luminosity 
function of Palomar 13 is adequately described by a rather 
flat mass function, with exponent x ~ 0.0 over the range 
21-5 ^ y ^ 24 (equivalent to 0.79 ^ M/Mq ^ 0.61; 
Bergbusch & VandenBerg 1992). If the initial mass func- 
tion of Palomar 13 had a more nearly universal form (for 
instance, near the "Salpeter value" oi x = 1.35; Salpeter 
1955), then the flat mass function wc find may be evi- 
dence for advanced dynamical evolution and the preferen- 
tial evaporation of low-mass stars. In short, the available 
data provide no reason to believe that the mass function of 
Palomar 13 is abnormal in any way, but firm conclusions 
must await deeper observations 

5.4. Non-Baryonic Dark Matter? 

In §5.2, it was found that, for an assumed mass-to-light 
ratio of Ty = 2 (i.e., a value appropriate for typical glob- 
ular clusters), the tidal radius of Palomar 13 predicted by 
equation (9) was in good agreement with the tidal radius 
of the low-concentration (c ^ 0.7) King-Michie model 
adopted by Siegel et al. (2001). As discussed in § 3.5 
and §4, it is possible to obtain an excellent match to the 
observed surface brightness and surface density profiles us- 
ing a King-Michie model which has c — 2.0 ± 0.15. In this 
case, the tidal radius of Palomar 13 is ~ 180 pc (with 
a formal uncertainty of ± 23% based on our Monte-Carlo 
simulations). From equation (9), we see that a tidal ra- 
dius of this size suggests a mass of M ~ 2 x 10^ Mq — 
i.e., considerably larger than our best-estimate for the dy- 
namical mass of Palomar 13 (assuming virial equilibrium) 
yet still far in excess of the expected mass. The actual un- 
certainties in the fitted tidal radii are difficult to gauge, 
but are unlikely to be smaller than ~ 50%, which trans- 
lates directly into a nearly a factor of two in mass. Thus, 
we conclude that the tidal radii inferred from the surface 
brightness and surface density profiles of Palomar 13 are: 
(1) highly dependent on how one chooses to model the 
profile; and (2) consistent with either the tidal disruption 
of an otherwise normal globular cluster, or the presence of 
substantial dark matter component in this object. 

In light of these results, it is natural to ask if the high 
velocity dispersion of Palomar 13 might also be the signa- 

^ Catalogs of pulsars and X-ray sources associated with globular 
Verbunt 2001). 



ture of a non-baryonic dark matter halo, similar to those 
belonging to Local Group dSph galaxies (Aaronson 1983; 
Mateo 1998; c./., Klessen k Kroupa 1998). This possibil- 
ity is particularly topical given recent claims from CDM 
simulations that the halos of large galaxies may contain 
large numbers of low-mass, "dark satellites" (e.g., Hi- 
rashita et al. 1999; Klypin et al. 1999; Bullock et al. 
2000). 

To examine this possibility, we consider a simple two- 
component model. First, we take the mass density profile 
of the luminous stellar component to be given by a Dehnen 
(1993) model, 

P*.™(r) = [(3-7)/4](Tp°,Aa3)(r/a)-'^(H-r/a)^-4 (10) 

where p'^ ; is the central luminosity density of this compo- 
nent, a is its scale radius, and T is its mass-to-light ratio. 
By contrast, the dark halo is assumed to have the form 
proposed by Navarro, Frenk & White (1997): 

p4r)=p"Ar/rs)-\l + r/r^)-' (11) 

where is the scale radius of the dark halo, and p° is the 
central dark matter density. Once these density distribu- 
tions are specified, the velocity dispersion profile follows 
from the the Jeans equation, 

dr r 

where (Trif) is the radial component of the velocity dis- 
persion and M{r) is the combined masses of the luminous 
and dark components. For our purposes, we assume that 
the velocity dispersion tensor is isotropic (i.e., /3 = 0). 
This theoretical velocity dispersion profile is then pro- 
jected onto the plane of the sky, and compared directly 
to the observed velocity dispersion. 

Before doing so, we require estimates of 7 and a from the 
surface brightness profile. The long-dashed curve in Fig- 
ure ^shows the best-fit Dchncn model, with 7 = 1.46±0.06 
and a = 135 ± 9" (15.9±1.0 pc), overlaid on the observed 
surface brightness profile. For an assumed mass-to-light 
ratio of Ty = 2, the total mass associated with this stel- 
lar component is M = 7 x 10'^ Mq. We also assume that, 
in projection, the dark matter matches that of the surface 
brightness profile, giving rg = 24 ± 1" (2.8±0.1 pc). This 
NFW model — fitted to the observed surface brightness 
profile — is shown by the solid curve in Figure M. The 
central mass density of this model is then left as a free pa- 
rameter which is adjusted until a reasonable match to the 
observed velocity dispersion is achieved. The correspond- 
ing velocity dispersion profile for a central mass density 
of p° = 80 ± 31 Mq pc~^ is shown as the solid curve in 
Figure 1^ (i.e., curve 3). This curve provides an adequate 
(though by no means unique) description of the observed 
velocity dispersion profile including the possible maximum 
at i? '--^ 1'. Although the mass of an NFW profile diverges 
as r — > cx), we note that over the range of our radial veloc- 
ity and surface brightness observations, the implied mass- 
to-light ratio of this two component model is Ty = 19^7 
— about half that found using a King-Michie model, but 

clusters list no such objects in Palomar 13 (Kulkarni & Anderson 1996; 
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still well above that expected for a normal stellar popula- 
tion. 

How does the central density of this putative dark halo 
compare to those found for the dSph satellites of the Milky 
Way? Table 1 of Mateo (1998) lists structural parame- 
ters for nine Galactic dSph galaxies having measured ve- 
locity dispersions: the central densities of these galaxies 
span the range 0.03 < p° < 0.6 Mq pc~'^, with a mean 
value of = 0.26 ± 0.07 Mq pc'^. Strictly speak- 
ing, these values refer to the combined densities of the 
dark and luminous components, p° = + How- 
ever, in view of the low stellar densities of these galaxies 
(e.g., p1„^ ^ 0.01 Mq pc~^), p° ~ p2 to good approxima- 
tion. Thus, the dark matter density needed to explain the 
measured velocity dispersion of Palomar 13 exceeds that of 
typical dSph galaxies by more than two orders of magni- 
tude. 

Following Kormendy (1990), we show in Figures |l^ and 
13 the scaling relations between central density, velocity 
dispersion, cr, absolute blue magnitude, Mb, core radius, 
Tc, and central _B-band surface brightness, /zb, for a sam- 
ple of Local Group dSph galaxies (Mateo 1998), dwarf 
irregular and spiral galaxies (Cote, Carignan & Freeman 
2000 and references therein) and Galactic globular clusters 
(Pryor & Meylan 1993; Harris 1996). Note that the densi- 
ties plotted here refer to those of the dark halos, with the 
exception of the globular clusters; for these objects, the 
densities are those calculated directly from the measured 
velocity dispersions. For Palomar 13, which is indicated by 
the large filled circle, we take the structural parameters for 
its putative dark matter halo to be those given by the best- 
fit NFW model. As these figures demonstrate, the central 
dark matter density in Palomar 13 is indeed much larger 
than those found for dwarf galaxies. On the other hand, it 
is also consistent with that expected from an extrapolation 
of galactic scaling relations to low luminosities. Indeed, it 
is remarkable that Palomar 13 falls so near the extrap- 
olations of the tight scaling relations for dwarf galaxies, 
although this agreement cannot be taken as unambiguous 
evidence for a common origin since Palomar 13 also falls 
near the edges of the more heterogeneous distributions for 
Galactic globular clusters (and, in the scaling relations in- 
volving fj,B, its location may be more consistent with the 
globular cluster distributions). Thus, the evidence from 
the observed scaling relations concerning the proper clas- 
sification of Palomar 13 as a bone fide globular cluster or 
an extreme dark-matter dwarf galaxy, is ambiguous. In 
any event, we note that a naive extrapolation of the dwarf 
galaxy scaling relations to globular cluster scales predicts 
central densities for dark matter halos in these objects that 
are sufficiently low (i.e., 1 <^ pj] ^ 100 Mq pc""^) to make 
their detection extremely challenging in all but the faintest 
clusters. In other words, for the vast majority of globular 
clusters, the inferred dark matter densities would fall well 
below the baryonic densities. 

Obviously, it would be of considerable interest to mea- 
sure velocity dispersions for additional faint, low surface 
brightness clusters. An especially interesting target is AM 
4, a cluster which shares many of the characteristics of 
Palomar 13: i.e., low luminosity (Ly ^ 400 Lv,q), low 
surface brightness {nv,o ~ 24.75 mag arcsec"^), and low 
metallicity ([Fe/H] ~ —2.0 dex). Like Palomar 13, it is 
also occupies a relatively isolated location in the Galactic 



halo, with Rq — 25.5 kpc. The central velocity dispersion 
predicted by equation (7) for a "normal" globular cluster 
of this surface brightness would be ap o <^ 0.3 km s~^, and 
it would be interesting to see if this expectation is borne 
out by radial velocity measurements. 

We consider three final scaling relations for dwarf galax- 
ies in Figure ^ In the upper panel of this figure, metal- 
licity is plotted against absolute magnitude for dE/dSph 
galaxies, along with the best-fit linear relation from Cote 
et al. (2000). Note that the relationship between metal- 
licity and luminosity constitutes a key difference between 
dwarf galaxies and globular clusters: i.e., galaxies span- 
ning a range of 10^ in luminosity follow a tight relation 
between metallicity and luminosity (e.g., Brodie & Huchra 
1991; Caldwell et al. 1992), whereas globular clusters obey 
no such correlation. While it is by no means clear that a 
simple extrapolation of this linear relation to My ^ —4 is 
justifiable, Palomar 13 is nevertheless considerably more 
metal-rich than might naively be expected for a dwarf of 
this luminosity. Furthermore, its CMD shows no evidence 
for an abundance spread, unlike Local Group dSph galax- 
ies which often show internal dispersions in metallicity of 

0.35 dex (e.g.. Cote et al. 2000; Shetrone, Cote & Sar- 
gent 2001). 

In the middle panel of Figure ^ mass-to-light ra- 
tio is plotted against absolute magnitude for all Local 
Group dE/dSph galaxies having measured velocity dis- 
persions. The dashed curve shows the relation expected 
for dwarf consisting of luminous stellar components with 
Tv = 2 that are embedded in dark matter halos of mass 
M = 2 X 10'^ Mq, as proposed by Mateo et al. (1993). 
The mass-to-light ratio predicted for Palomar 13 by this 
relation would he Ty 7000, making it an extreme out- 
lier. On the other hand, Palomar 13 does seem to fol- 
low the correlation between mass-to-light ratio, metallic- 
ity and central surface brightness for Local Group dwarfs, 
shown in the lower panel of Figure |l^ (Prada & Burkert 
2001). We can only conclude that the evidence regarding a 
possible connection between Palomar 13 and dark-matter- 
dominated dwarf galaxies is ambiguous and that, in light 
of the complexity of gas cooling, star formation and feed- 
back in low-mass systems (see, e.g.. Figure 8 of Mateo 
1998), a definitive conclusion is not possible at the present 
time. 

5.5. Modified Newtonian Dynamics? 

A much different interpetation of the large velocity dis- 
persion of dSph galaxies has been proposed by Milgrom 
(1983a), who suggested that at low accelerations, a ^ oq, 
Newtonian gravity should be revised to include a repulsive 
term. If this modified theory of Newtonian dynamics (i.e., 
MOND) is correct, then this there would be no need to 
invoke dark matter as an explanation for the high veloc- 
ity dispersions of dwarf galaxies (Milgrom 1983b; 1995). 
But if MOND is to be a viable alternative to dark mat- 
ter, then it must apply to all systems in which the the 
acceleration falls below the characteristic MOND value, 
qq ~ 1.2 X 10~* cm s~^ — value chosen to match the ro- 
tation curves of low-surface brightness galaxies (Begeman, 
Broeils & Sanders 1991; McGaugh & de Blok 1998). Since 
this condition is met for Palomar 13, we now examine its 
internal dynamics within the MOND framework. 

There are two relevant accelerations for Palomar 13: the 
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acceleration arising from the internal mass distribution, 
and an external component due to the Milky Way. We ap- 
proximate the internal acceleration as Oint = GM{r)/r'^, 
where M(r) is mass enclosed with a radius, r, for the 
best-fit King model shown in Figure 0. For an assumed 
mass-to-light ratio of Ty = 2, the internal acceleration 
reaches a maximum value of ai„t ~ 1 x 10^^ cm s^^ at 
r ~ Tc- The external acceleration may be approximated 
as ttext = QI/Rgi or aoxt = 6 x 10~^ cm s~^ for 6c = 220 
km s-i and Rq = 25.3 kpc. Thus, according to MOND, 
Palomar 13 is in the "quasi-Newtonian weak- field limit" 
since a-mt < aoxt < ao- In this regime, the Newtonian 
mass is related to the MOND mass by a factor, fj,, which 
falls in the range 1 ^ ^ ao/ooxt (Gerhard & Spcrgcl 
1992). From Table 4, the best-fit King-Michie model has 
Tv — 40l^y in solar units. Thus, the mass-to-light ratio 
in MOND would be reduced by — at most — a factor of 
ao/flcxt, giving Ty ~ 20^9^. Apparently, the high mass- 
to-light ratio of Palomar 13 persists (at the 2a level) even 
within the MOND framework. We conclude that this ex- 
planation appears less appealing than those involving tidal 
disruption or dark matter. 

6. SUMMARY AND CONCLUSIONS 

An imaging and spectroscopic survey of globular clus- 
ters in the outer halo of the Milky Way has shown Palomar 
13 to exhibit several puzzling features. 

From an abundance analysis of our HIRES spectra for 
a single RGB star and isochrone fitting to the CMD, 
we estimate the metallicity of Palomar 13 to be [Fe/H] 
= — 1.9 ± 0.2 dex. We find a true distance modulus of 
(m — Af)o = 16.93 ± 0.1 mag and a distance of Rs — 
24:.3^\'i kpc, placing Palomar 13 squarely in the outer 
halo. With My = —3.8 mag, this is the lowest luminosity 
object in our survey. The sparsely populated red giant and 
horizontal branches preclude a precise age determination, 
but a 14 Gyr isochrone is found to provide a reasonable 
match to the observed fiducial sequences. Thus, in terms 
of its stellar populations, Palomar 13 appears typical of 
globular clusters in the outer halo. 

However, from a sample of 21 probable members with 
radial velocity measurements from HIRES, we find the sys- 
temic velocity and intrinsic velocity dispersion of Palomar 
13 to be {vr)s = 24.1 ± 0.5 km and CTp = 2.2 ± 0.4 
km s^^, respectively. This velocity dispersion is roughly 
four times higher than that expected on the basis of the 
observed luminosity and concentration. Taken at face 



value, this dispersion suggests a mass-to-light ratio of 
Ty — 40t^7, '^ith the exact value depending on the de- 
tails of the modeling. This mass-to-light ratio is far higher 
than those found for globular clusters, but is comparable to 
those of Galactic dSph galaxies. The surface density pro- 
file is also found to be anomalous among Galactic globular 
clusters, in that it shows a significant number of apparent 
member stars beyond the tidal radius of a low concentra- 
tion King-Michie model (i.e., c ~ 0.7, as has usually been 
claimed for Palomar 13). On the other hand, an excellent 
fit to the surface brightness profile over all radii is provided 
by a King-Michie model with c ^ 2, as well as by models 
of the type suggested by Dehnen (1993) and Navarro et al. 
(1997).' 

Based on the measured velocity dispersion and surface 
density profile, two possibilities present themselves: (1) ei- 
ther Palomar 13 has recently undergone catastrophic heat- 
ing and mass loss, as might be expected on the basis of its 
orbital properties; or (2) its spatial extent, central con- 
centration, and mass have been greatly underestimated in 
previous studies. Other factors, such as a large fraction 
of spectroscopic binary stars or the presence of a large 
number of heavy remnants, may play contribute to Palo- 
mar 13's unusual properties but are, on their own, unlikely 
to provide a complete explanation of the observed struc- 
tural and kinematic properties. Additional observations 
— particularly deep imaging with HST and ground-based 
mosaic cameras, and radial velocity measurements for an 
expanded sample of stars at large radii — should help dis- 
criminate between the various possibilities. In any event, 
it is clear that Palomar 13 merits further study, as it of- 
fers new insights into the structure and formation of the 
Galactic halo. 
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Table 1 
Observing Log 



Telescope Date Spectrograph Imager 



Keck II 10/09/1999 ... LRIS + SF/ filters 

CFHT 14/07/1999 ... CFH12K + VR filters 

Keck I 21-22/08/1998 HIRES + CI Decker 

Keck I 15-16/10/1998 HIRES + C5 Decker 

Keck I 15-16/07/1999 HIRES -t- CI Decker 

Keck I 11/08/1999 HIRES + CI Decker 
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Table 2 

Radial Velocities for Candidate Red Giants in Palomar 13 



ID 


R 

(") 


V 
(mag) 


[B-V) 
(mag) 


T 
(sec) 


HJD 

2,450,000+ 


Rtd 


(k: 


m s^^ 


) 


(k: 


{Vr) 

m s "'^ 


) 


p(m) 

(%) 


Memb 


ORS- 


13 


113 


16.55 


0, 


.89 


180 


1046.9063 


17.61 


-38, 


.38±0, 


.48 


-38. 


.90±0. 


.29 





N 














240 


1047.9159 


24.72 


-39, 


.17±0, 


.35 












ORS- 


■118 


102 


17.00 


0, 


,93 


270 


1046.9147 


13.04 


25, 


.43±0, 


.64 


24. 


.92±0. 


.21 





Y 














600 


1374.9741 


27.43 


25, 


.07±0, 


.32 
























720 


1375.9423 


29.27 


24, 


.68±0, 


.30 












ORS- 


■6 


125 


17.03 


1, 


.10 


240 


1046.9104 


18.00 


0, 


.68±0. 


.47 


0. 


.68±0. 


.47 





N 


ORS- 


■103 


134 


17.04 


1, 


.00 


360 


1046.9196 


24.48 


-56, 


.72±0, 


.35 


-56. 


.72±0. 


.35 





N 


ORS- 


■110 


100 


17.21 


1, 


.03 


600 


1046.9266 


23.96 


5, 


.40±0, 


.36 


5. 


.40±0. 


.36 





N 


ORS- 


■72 


60 


17.64 


0, 


.85 


240 


1101.7681 


10.22 


28, 


.64±0, 


.80 


28. 


.79±0. 


.27 


88 


Y 














750 


1046.9350 


13.45 


29, 


.30±0, 


.62 
























750 


1047.9235 


13.14 


30, 


.05±0, 


.64 
























900 


1374.9973 


24.38 


28, 


.26±0, 


.36 












ORS- 


■31 


95 


17.76 


0, 


.83 


750 


1046.9456 


10.70 


26, 


.31±0, 


.77 


25. 


.09±0. 


.35 


79 


Y 














900 


1374.9844 


22.02 


24, 


.78±0, 


.39 












ORS- 


■91 


27 


17.81 


0, 


.67 


750 


1046.9554 


6.72 


24, 


.46±1, 


.17 


24. 


.59±0. 


.60 


99 


Y 














900 


1375.0091 


12.39 


24, 


.64±0, 


.67 












V2 







17.91 


0, 


.49 


1200 


1375.1025 


9.32 


25, 


.39±0, 


.87 


25. 


.39±0. 


.87 




Y 


ORS- 


■32 


62 


18.02 


0, 


.89 


900 


1046.9664 


11.56 


20, 


.46±0, 


.72 


19. 


.68±0. 


.26 


35 


Y 














900 


1375.1170 


18.95 


19, 


.31±0, 


.45 
























1200 


1376.1122 


18.60 


18, 


.69±0, 


.46 
























1200 


1401.9353 


17.23 


20, 


.90±0, 


.49 












ORS- 


■41 


64 


18.59 


0, 


.76 


320 


1101.7397 


3.91 


18, 


.37±1 


.84 


19. 


.24±0. 


.47 


88 


Y 














1200 


1046.9987 


8.72 


18, 


.91±0, 


.93 
























1500 


1375.9567 


15.56 


19, 


.43±0, 


.54 












ORS- 


■1 


151 


18.59 


0, 


.82 


1200 


1046.9799 


13.35 


-37, 


.48±0, 


.63 


-37. 


.48±0. 


.63 





N 


ORS- 


■101 


39 


18.65 


0, 


.94 


1200 


1047.0141 


18.88 


-17, 


.39±0, 


.45 


-17. 


.39±0. 


.45 





N 


ORS- 


■14 


92 


18.76 


0, 


.74 


1500 


1047.0337 


17.63 


-55, 


.31±0, 


.48 


-55, 


.31±0, 


.48 





N 


ORS- 


■86 


14 


18.81 


0, 


.77 


1800 


1047.8735 


10.01 


24, 


.47±0, 


.82 


24, 


.47±0, 


.82 


99 


Y 


ORS- 


■36 


41 


18.98 


0, 


.75 


1800 


1047.0760 


9.13 


25, 


.29±0, 


.89 


25, 


.29±0. 


.89 


93 


Y 


ORS- 


■63 


6 


19.02 


0, 


.76 


1800 


1047.0530 


9.12 


26, 


.09±0, 


.89 


25, 


.09±0. 


.42 


62 


Y 














2100 


1375.0805 


18.35 


24, 


.82±0, 


.47 












ORS- 


■87 


16 


19.05 


0, 


.72 


1800 


1047.9021 


9.06 


26 


.55±0, 


.90 


26, 


.55±0, 


.90 


94 


Y 


910 




111 


19.28 


0, 


.73 


600 


1101.7614 


4.15 


23 


.77±1, 


.75 


23, 


.77±1, 


.75 


89 


Y 


911 




52 


19.36 


0, 


.75 


600 


1101.7306 


3.86 


23, 


.19±1, 


.86 


23, 


.19±1. 


.86 


99 


Y 


ORS- 


■23 


77 


19.49 


1, 


.50 


2400 


1047.1068 


10.01 


18, 


.57±0, 


.82 


18, 


.57±0. 


.82 





N 


915 




9 


19.58 


0, 


.73 


1800 


1375.0284 


11.68 


24, 


.69±0, 


.71 


24, 


.69±0, 


.71 




Y 


ORS- 


■18 


112 


19.62 


0. 


.73 


1000 


1102.7253 


4.36 


25, 


.16±1, 


.68 


25, 


.16±1. 


.68 


98 


Y 


ORS- 


■38 


40 


19.65 


0, 


.74 


1000 


1101.7489 


4.49 


19, 


.52±1, 


.64 


21. 


.64±0. 


.72 


88 


Y 














2400 


1376.0856 


10.57 


22 


.11±0, 


.78 












931 




18 


19.70 


0, 


.72 


2100 


1375.0528 


11.70 


25, 


.45±0, 


.71 


25, 


.45±0. 


.71 


91 


Y 


ORS- 


■78 


19 


19.73 


0, 


.72 


1000 


1102.7404 


5.26 


22, 


.34±1, 


.44 


22, 


.34±1. 


.44 


78 


Y 


ORS- 


■50 


22 


19.76 


0, 


.72 


2700 


1376.0226 


9.66 


24, 


.15±0, 


.85 


24, 


.15±0. 


.85 


98 


Y 


ORS- 


■96 


32 


19.80 


0, 


.71 


1000 


1102.7542 


6.52 


25, 


.38±1, 


.20 


25, 


.38±1. 


.20 


99 


Y 


ORS- 


■5 


117 


19.86 


0, 


.72 


2700 


1375.9867 


9.99 


20, 


.95±0, 


.82 


20, 


.95±0, 


.82 


98 


Y 


ORS- 


■88 


37 


20.11 


0, 


.71 


2100 


1376.0541 


7.34 


24, 


.61±1, 


.08 


24, 


.61±1. 


.08 


78 


Y 
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Table 3 

Observed and Derived Parameters for Palomar 13 



Parameter 


Symbol 


Value 


Units 


Reference^ 


Right Ascension 


a(J2000) 


23:06:44.8 


h:m:s 


1 


Declination 


(5(J2000) 


-1-12:46:18 


o:/:// 


1 


Apparent Distance Modulus 


(m - M)v 


17.27±0.1 


mag 


2 


True Distance Modulus 


(m ^ M)o 


16.93±0.1 


mag 


o 
z 


Reddening 


E{B - V) 


0.11±0.02 


mag 


3 


Galactic Coordinates 


il,b) 


(87.1,-42.7) 


deg 


-1 
1 


Distance 


Rs 


24.3tJ;i 


kpc 


2 


Galactocentric Distance 


Rg 


25.3tJ-i 


kpc 


2 


Metallicity 


[Fe/H] 


-1.9 ±0.2 


dex 


2 


Perigalacticon distance 


Rp 


11.2 


kpc 


1 


Apogalacticon distance 


Ra 


80.2 


kpc 


1 


Orbital Period 


P 


1.1 


Gyr 


1 


Eccentricity 


e 


0.76 




1 


Maximum Likelihood Mean Velocity 


{Vr)s 


24.1 ±0.5 


km s 


2 


Maximum Likelihood Dispersion 


a„ 


2.2 ±0.4 


km s 


2 


King-Michie Model of Siegel et al. (2001 ) 






Core Radius 


rc 


39 


arcsec 


1 






4.6 


pc 


1 


Tidal Radius 


n 


188 


arcsec 


1 






23 


pc 


1 


Central Projected Luminosity Density 




23.39 


mag arcsec 


rt 

2 






15 


Lv^Q pc 


r> 

z 


Concentration 


c 


0.7 




-1 

1 


Integrated Luminosity 


Lv 


1.2 X 10^ 


Lv,Q 


1 


Best-fit King-Michie Model 






Core Radius 




14±2 


arcsec 


2 






1.70±0.24 


pc 


2 


Tidal Radius 


n 


26±6 


arcmin 


2 






182±41 


pc 


2 


Central Projected Luminosity Density 




22.54t°-20 


mag arcsec"^ 


2 






33.8±5.8 


Lv.Q pc 


2 


Concentration 


c 


2.0±0.15 




2 


Integrated Luminosity 


Lv 


(2.8 ±0.4) X 10^ 




2 


Mass-to-Light Ratio 


Tv 


40124 


Mq/Lvq 


2 




Dehnen Model 






Power-Law Index 


7 


1.46±0.06 




2 


Scale Radius 


a 


135±9 


arcsec 


2 






15.9±1.0 


pc 


2 


Integrated Luminosity 


Lv 


(3.5 ±0.1) X 10^ 


Lv,Q 


2 




NFW Model 






Scale Radius 




24±1 


arcsec 


2 






2.8±0.1 


pc 


2 


Central Luminosity Density 




4.3±0.4 


LvQ pc"^ 


2 


Central Mass Density 




80±31 


Mq pC-3 


2 


Central Mass-to-Light Ratio 




191? 


Mq/Lvq 


2 



^References for Table | (1) Siegel et al. (2001); (2) This paper; (3) Schlegel, Finkbeiner & Davis 
(1998) 
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Fig. 1. — V^-band image of Palomar 13 taken with the Low-Resolution Imaging Spectrometer on the Keck II telescope. This image measures 
5f8x7f3. The solid circle shows the core radius (rc = 14") of the best-fit King-Michie model. The dashed circle denotes the "break" radius 
where the possible excess in the background-subtracted surface density profile begins: Ri, ~ 1' (see §5.1). The shortest of the three arrows 
shows the direction of the absolute proper motion from Siegel et al. (2001). The intermediate arrow points in the direction perpendicular to 
the Galactic plane, while the long arrow points in the direction of the Galactic center. Stars judged to be probable cluster members based 
on their radial velocities and proper motions are circled. The size of each circle is proportional to the absolute value of the difference in 
the measured radial velocity from the systemic velocity of the cluster. Red and blue circles denote stars with positive and negative velocity 
residuals, respectively. Non-members are indicated by open squares. 
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Fig. 2.- (Upper Panels) Difference in V magnitudes for stars in common between our LRIS photometric study and those of Ortolani, 
Rosino & Sandagc (1985; 085), Borissova, Markov & Spassova (1997; B97) and Siegel ct al. (2001; SOI). (Lower Panels) Difference in (B-V) 
color for stars in common between our LRIS photometric study and those of Ortolani, Rosino & Sandage (1985; 085), Borissova, Markov & 
Spassova (1997; B97) and Siegel et al. (2001; SOI). 
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Fig. 3. — Proper motion membership probabiUty from Siegcl ct al. (2001), plotted against mean radial velocity for candidate Palomar 13 
stars. The dotted curves show approximate velocity distributions along this line of sight for disk and halo field stars (see text for details). 
The dashed vertical line indicates the systemic velocity of Palomar 13. The error bar in the upper left corner showrs the typical uncertainty 
in the measured radial velocity. Open squares show non-members, while probable cluster members are indicated by the filled circles. Note 
that this latter sample includes one star (ORS-118) which has P(/i) = 0% according to Siegel et al. (2001). We argue that this star is in fact 
a cluster member since it: (1) has a radial velocity which is indistinguishable from cluster systemic velocity; (2) is located precisely on the 
red giant branch in the cluster color-magnitude diagram; and (3) has a metallicity of [Fe/H] = —1.98 ± 0.31 dex measured from our HIRES 
spectra. 
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Fig. 4. — (Left Panel) V, {B — V) color magnitude diagram for Palomar 13 based on our LRIS images. A total of 840 objects are plotted 
in this figure. Open circles indicate probable cluster members, determined from our HIRES radial velocities. The size of each circle is 
proportional to the absolute value of velocity residual with respect to the cluster's systemic velocity. Radial velocity non-member stars are 
indicated by the filled squares. Known RR Lyra« variables are indicated by the stars; the circled variable is V2. The solid curve shows an 
isochrone from Bergbusch & VandenBerg (1992) having an age of T = 14 Gyr and a metallicity [Fe/H] = —1.78, which has been shifted by 
E{B -V) = 0.11 mag (Schlegel et al. 1998) and (m-M)v = 17.27 mag. (Right Panel) V, (V - I) color-magnitude diagram for Palomar 13 
derived from our LRIS images. The symbols are the same as in the previous panel. 
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Fig. 5. — Surface density profile for Palomar 1,3 based on our Kcclc and CFHT photometry (squares and circles, respectively). Open 
symbols show the surface density profile before background subtraction; filled symbols show the profile found after subtracting backgrounds 
of Sf, = 0.65 and 0.62 stars per arcmin^^, respectively, for the Keck and CFHT observations. For comparison, the surface density profile of 
Siegel et al. (2001) is illustrated by the open triangles, after scaling their counts upward by a factor of (629/119) ~ 5.28 (i.e., the ratio of the 
rmmber of stars used to derive the respective profiles). The dotted curve shows a King model with concentration parameter c = 0.7 and core 
radius r'c = 39" suggested by Siegel et al. (2001). The dashed curve shows the formal best-fit King-Michie model. The three dashed lines 
beginning at i? ~ 1' indicate profiles of the form predicted by Johnston et al (2001) for extra-tidal stars: S oc R ' with 7 = 0.5, 1.5 and 2.5. 
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Fig. 6. — Instrumental V, (V — r) color magnitude diagrams based on our CFHT images in three radial regimes. The left panel shows 
all objects within the canonical tidal radius, rt = 188", reported by Siegel et al. (2001). The dotted rectangles shows the magnitude and 
color limits used to reject probable non-members. Circled points indicate probable radial velocity members for which we have reliable CFHT 
photometry. All objects located beyond this radius are plotted in the middle panel, where the dashed curve shows the Palomar 13 fiducial 
sequence from the previous panel. Objects at radii greater than twice the canonical tidal radius are plotted in the right panel. 
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Fig. 7. — V-band surface brightness profile for Palomar 13. Pentagons and squares represent surface photometry and star counts based 
on our Keck imaging, respectively. Surface photometry and star counts using our CFHT imaging are indicated by the triangles and circles. 
The King-Michie model which best fits this surface brightness profile is shown by the dashed curve. By contrast, the King-Michie model with 
structural parameters found by Siegel et al. (2001) from their surface density profile (see Figure 5) is indicated by dotted curve, shifted in 
the vertical direction. The solid curve shows the NFW model which best fits the surface density profile, while the long-dashed curve shows 
the best-fit Dehnen model. The dotted curves at the cluster center shows typical stellar profiles on our Keck and CFHT images, scaled to the 
luminosity of V2 — the bright variable which defines the cluster center. 
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Fig. 8. — (Upper Panel) Heliocentric radial velocity for probable cluster members plotted against distance from the cluster center (filled 
circles). The dashed line indicated the cluster's systemic velocity. The arrows denote the approximate "break" radius, Rf, ~ 1', where the 
possible excess in the cluster's surface density profile begins. (Lower Panel) Velocity dispersion profile for Palomar 1,3. The open circle shows 
the intrinsic velocity dispersion for the entire sample of stars, plotted at the mean radius of the sample. The two filled circles show the 
intrinsic velocity dispersion which is found if the sample is divided into two radial bins containing 11 and 10 stars each. The diflferent curves 
(labeled 1-4) indicate various model fits to the data; a detailed description of each model is presented in the text. 
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Fig. 9. — (Upper Panel) Histogram of position angles for stars in our photometric database which are located in the radial range 1' < R < 
2'.88. The dashed line shows the mean number of stars per azimuthal bin, while the dotted curve shows the best-fit sinusoid: A'^ oc (4.5 ± 
1.7) sin [20 + (60° ± 11°)]. This curve has density maxima at &p = 15° and 195°. (Lower Panel) Mean stellar surface density in the range 
1' < /? < 2.'88 for each of four 90° sectors centered on cluster center. The dashed and dotted curves correspond to the mean surface density 
and sinusoid described in the upper panel (i.e., with a best- fit amplitude of 2.4it0.9 arcmin"^). Surface densities are shown for three different 
orientations of the quadrants: in the first case (filled circles), the sector boundaries are offset by ±45° from the position angle of the vector 
which points in the direction of Palomar 13's absolute proper motion (i.e., the short arrow in Figure hh. In the second case (open circles), the 
boundaries^ are offset by ±45° from the position angle of the vector which points in the direction of the Galactic center (i.e., the long arrow 
in Figure hi). In the final case (open squares), the boundaries are offset by ±45° fmm the position angle of the vector which points in the 
direction perpendicular to the Galactic disk (i.e., the intermediate arrow in Figure hi). The orientations of these three vectors are indicated 
by the respective arrows. There is only marginal evidence for a preferred orientation of Palomar 13 stars in this radial range. 
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Fig. 10. — (Upper Panel) Mean radial velocity vs. position angle, 0, for the sample of 21 probable member stars. The dashed line indicates 
the systemic velocity, while the dotted curve shows the best-fit sinusoid of the form (vr) = (vr)s + Asin{0 + 0p), where {vr)s = 24.1 km . 
A weighted least-squares fit yields A = 2.1 ± 1.0 km and 0p = 48 ± 47°, while an unweighted fit gives A = 0.8 ± 1.0 km s"! and 
0p = 104 ± 130° . We conclude that there is no evidence for rotation among the full sample of Palomar 13 members. (Lower Panel) Same as 
above, except for the eight stars located at, or beyond, the "break radius" of Rf, ~ 1'. The weighted best-fit sinusoid is shown by the dotted 
curve, which has parameters A = 3.0 ± 3.0 km and 0p = 60 ± 65°. An unweighted fit gives A = 2.0 ± 1.7 km and 0p = 84 ± 44°. 
Thus, as with the full sample, we find no evidence for statistically significant rotation for those stars at large radii. 
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Fig. 11. — Comparison of the observed luminosity function for Palomax 13 with theoretical luminosity functions for V ^ 21.5 mag. 
Data points from Keck are shown by the squares; circle indicate measurements from CFHT. The theoretical luminosity functions have been 
generated from the [Fe/H] = —1.78 dex, 14 Gyr isochrone of Bergbusch & VandenBerg (1992). Six different values of the mass function 
exponent, x, are shown; from top to bottom, the dashed curves have x = 2.5, 2.0, 1.5, 1.0, 0.5 and 0.0. The steep dropoffs at V ^ 23.5 mag 
(CFHT) and F ^ 24 mag (Keck) are due to photometric incompleteness. Thus, to the completeness limits of our data, we find no evidence 
for a steeply rising luminosity function in Palomar 13. On the contrary, the observed luminosity function appears quite fiat, with x ~ 0, 
indicating an possible underabundance of low-mass stars. 
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Fig. 12. — Scaling relations for stellar systems. Local Group dwarf spheroidal galaxies (data from Mateo 1998) are indicated by open circles. 
Open squares show the results for a sample of dwarf irregulars and spirals taken from Cote et al. (2000). Galactic globular clusters having 
measured velocity dispersions are shown by the dots (Pryor & Meylan 1993; Harris 1996) . Palomar 13 is indicated by the large filled circle. 
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Fig. 13. — Additional scaling relations for galaxies and Galactic globular clusters. The symbols are the same as in the previous figure. 
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Fig. 14. — (Upper Panel) Metallicity-luminosity relation for dE/dSph galaxies (filled circles). Palomaj: 13 is indicated by the open square. 
The dashed line shows the best- fit relation from Cote et al. (2000), extrapolated to faint magnitudes. (Middle Panel) V-band mass-to-light 
ratio, T\A, versus absolute visual magnitude for Local Group dE/dSph galaxies (see Cote et al. 1999). The dashed curve shows the relation 
for dwarf consisting of luminous stellar components with = 2 that are embedded in dark matter halos of mass M = 2 x 10^ A^q. Palomar 
13 obeys neither this, nor the above, empirical scaling relation for dwarf galaxies. (Lower Panel) Correlation between mass-to-light ratio, 
metallicity and central V-band surface brightness for Local Group dE/dSph galaxies (filled circles). The dashed line shows the linear relation 
proposed by Prada & Burkert (2001). Palomar 13 is in reasonable agreement with the proposed correlation. 



This figure "fig01.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/0203410vl 



